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Abstract 

, The response of a massive body to gravitational waves is described on the 

■ microscopic level. The results shed a new light on the commonly used oscil- 

lator model. It is shown that apart from the non-resonant tidal motion the 
energy transfer from a gravitational wave to an electromagnetically coupled 
body is in general restricted to the surface, whereas gravitational coupling 
\ gives rise to bulk excitation of quadrupole modes, but several orders of mag- 

. nitude smaller. A microscopic detector making use of the effect is suggested. 

CO 

O ■ 1 Introduction 

Gravitational waves were already considered by Einstein as the wave solutions of 
the linearized field equations of gravity. There is indirect evidence of their exis- 
tence through systems of binary pulsars that loose energy in form of gravitational 
radiation [jjj, their direct experimental measurement presently is one of the most 
challenging tasks in gravitational physics. Very sensitive detectors operating at the 
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£~ ' ■ quantum limit are needed to detect directly gravitational waves from cosmic events 

such as collapsing or colliding star systems. There are basically two different types of 
detectors: resonant mass antennas based on the resonant excitation of quadrupolc- 
k> \ type modes of a appropriately chosen massive body, like the bar detectors conceived 

j_j ■ by Weber ||] , and laser interferometric devices that detect the direction-dependent 

variation of the proper distance between the mirrors of a Michelson interferometer 
H . Detectors of both types are presently under construction j| . 

Commonly a resonant mass antenna is described in Riemannian normal coordi- 
nates with respect to its center of mass, the proper frame of reference (PFR). The 
detector is analyzed in term of normal modes, idealized by a spring that couples 
two masses Q , the resonant energy input is calculated. The intention of the work 
presented in this article was twofold: First, to validate the results obtained from 
the normal mode model by microscopic considerations, second to give a comple- 
mentary description of the detector in the reference system of the wave, in which 
the linearized solutions of the Einstein field equations are computed. Our model 
is based on the local properties of a detector in terms of the fundamental binding 
forces, electromagnetic and other, which we consider in both the PFR and the wave 
system. Analyzing the deviations from the tidal motion, we find that the energy 
input from gravitational waves on an electromagnetically coupled massive body is 
restricted to the surface of the body, whereas gravitational coupling leads to true 
bulk excitation of quadrupole modes. This result does not contradict the normal 
mode picture at all, rather it presents a complementary viewpoint that has eluded 
the normal mode analysis. The reason is that though the energy input into any 
single mode is nonlocal, in the special case of gravitational waves the superposition 
of all excited modes describes a localized excitation. Based on our observations we 
propose a new, microscopic type of detector. 
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2 Frames of reference 



There exist two frames of reference that can be used for the analysis of gravitational 
waves and gravitational wave detectors. The natural system to study the waves is 
a perturbed Minkowski system for which the linearized Einstein field equations 
are solved. In this reference system the plane wave is gauged, conventionally a 
transverse-traceless (TT) gauge is chosen [g . On the other hand, the use of the PFR 
system with Riemannian normal coordinates with respect to the center of mass is the 
natural system for the study of a detector. For our purpose to study of the detector 
on the microscopic level, both system have advantages and disadvantages. In the 
PFR system we have common argument || that in the metric of the gravitational 
wave field leads to variations of the electromagnetic field of order 

SA/A ~ {L 2 /X 2 ) h TT (1) 

where L is the distance from the origin, A and h TT the wave length and amplitude 
of the gravitational wave, respectively, whereas the tidal gravitational forces lead to 
displacements of the constituent particles of the body that in turn lead to variations 
of the electromagnetic field of order 

SA/A ~ h TT (2) 

so that the metric effect can be neglected. This means that in the PFR system the 
unperturbed solutions of the Maxwell equations can be used. On the other hand, in 
the PFR system we have to deal with the general problem of general relativity that 
local energy densities have no invariant meaning, so that we cannot easily control 
the exchange of energy between the wave and the detector. Here the TT system 
has the advantage that we can employ a Hamiltonian description of the constituent 
particles, and because the system is in principle a special relativistic one, we have 
the standard laws of energy-momentum conservation. Therefore the TT system 
proves useful for our considerations. Because considerations in both coordinate 
system contribute to the understanding of the detector response, we formulate our 
results in both systems. 

We first show that in the TT system no energy is transferred to a system of 
non-interacting particles; in the PFR system this corresponds to the well-known 
quadrupole type oscillations. Next we extend the discussion to electromagnetically 
coupled systems, where we can use the standard Coulomb force in the PFR system, 
whereas we have to solve the metric Maxwell equations in the TT system. Finally 
we discuss a gravitationally coupled system and show that there exists a decisive 
difference in the response of the detector which is due to the different nature of the 
fundamental forces. 



3 Particle Motion 

A point-like test mass with electric charge e can be described in the presence of 
electromagnetic fields and gravitation by the Hamiltonian JtJ 

H = c^[p tl -eA^)g^{p l/ -eA v ) (3) 

where we use coordinates x^, ji = 0, 1, 2, 3, 4 for space-time with metric g^ix*") and 

signature H , p^ are the momentum coordinates in the cotangent space, and 

A p (a; K ) is the electromagnetic four-potential. We also use the notation (ct,x,y,z) 
in an obvious manner. The evolution parameter will be denoted by r. The Hamilto- 
nian is conserved, dH/dr = 0, it represents the rest mass m — H/c 2 oi the particle. 
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The canonical equations of motion are given by 



dH 



dH 



(4) 



The constancy of H is equivalent to x >1 g >lv x v 
evolution parameter is the proper time. 

A gravitational wave propagating in z-direction with 
modes is described in TT gauge by the metric tensor 
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where h^v is a small perturbation of the Minkowski metric r)^ u @J . This perturbation 
acts as classical, special relativistic field, so that the system can be treated within 
special relativity, except that the interpretation of the energy must be treated with 
care [[|. In the absence of an electromagnetic field we obtain for this metric a 
Hamiltonian that leads to four conserved quantities: H,p x ,p y , and po+p 3 = E/c + 
p z where E is the energy of the particle in the sense of special relativity. Note 
that p 3 = —mi, so that the difference between the energy and the conventional 
z-momentum is conserved. This is natural, since the gravitational wave not only 
carries energy E w , but also momentum P w with the relation E w — c \P W \ that holds 
for all massless objects in special relativity. The exchange AE of energy between 
the wave and a test mass thus is always accompanied with an exchange AP of 
momentum: 



AE = cAP. 



(6) 



The existence of four conserved quantities now allows us to integrate the equations 
of motion completely: 



• n l • i 

x = tPo x 

xl = h. (d 1 1 Pi + 9 12 P2) i 3 



±(9 21 Pi+3 22 P2) 



pi = const 

1 , ,1m 2 

Pa = o Po+P3j + ^7 — - — r 

2 2(po+p 3 ) 

1 , ,1m 2 

P3 = o P0+P3 -^7 ; r 

2 2(po+p 3 ) 



P2 = const 

\pa9 ab Pb 



(7) 



c 2 - 



rn 
1 

rn 



2Pa9 Pb 



(8) 



where we use the indices a, b for a summation over 1, 2 only. Let us consider a wave 
pulse. We denote the initial conditions before the arrival of the pulse by and 
have 
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where h ab = g ab — r\ ab . Thus after the pulse, where the perturbation h is zero 
again, the particle has the same four-momentum as before, and the only possible 
effect is a displacement of the straight trajectory after the pulse from the one before 
the pulse. If a particle is initially at rest, it stays at rest in this reference frame. 
Two particles that are at rest relative to each other, remain at rest relative to each 
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other, though the proper distance between them changes with the wave amplitude. 
Thus free test particles do not take up energy from a gravitational wave. When we 
restrict this statement to a comparison of the energy before and after a wave pulse 
or wave train, we are moreover free from the ambiguity of the energy definition in 
general relativity. 

We now look at this result in the PFR system. In the PFR system particles 
do not stay at rest in the wave field, but move under tidal accelerations. In PFR 
coordinates (t,x,y,z) the tidal accelerations produced by the wave are 

^ a 6 ~b ld 2 h ab „ b 

RaObQX = 7;—^T x ( 10 ) 



dP " uu " 2 dP 

For a particle that is initially at rest at x a ^ before a wave pulse arrives, the solution 
is simply given by 

r(i)= x ^ + h ab (t)x b ^. (ii) 

(Note that R a obo and R a obo differ only by terms of order O (h 2 )). This is the well- 
known quadrupole-like tidal motion of a system of particles around the origin of 
the frame of reference. The coordinate distance of the particles varies according 
to the changing metric distance between them. We can ascribe a standard kinetic 
energy to this motion, but there will be additional contributions to the conserved 
energy from the metric. Our analysis in the TT system has shown that the energy 
of the wave does not vary in this case. The coordinate transformation from the TT 
to the PFR system is time-dependent, so that the energy conservation in the TT 
system translates into a more complicated conservation law in the PFR system. 

Because free particles do not effectively take up energy from a gravitational wave, 
we have to take the coupling between particles into account in order to describe the 
response of a detector. The coupling is basically of electromagnetic nature in small 
massive bodies, but may also be of gravitational nature in large bodies. Because 
the effects of the gravitational wave are so small, massive detectors must be cooled 
to zero temperature as near as possible. Thus the ground state of a body where 
the constituent atoms or ions are at rest relative to each other except for quantum 
effects will serve as an appropriate model. Rotational motion with respect to the 
TT frame of reference must be taken into account. Internal motion, thermal or 
otherwise, of the atoms leads to forces of order mAvdh/dt; we do not consider 
this kind of phonon-graviton interaction in this article, because it is the idea of the 
Weber detector to excite the quadrupole modes, not to enhance already excited 
modes. 



4 Electromagnetic Field 

We first look at the Coulomb potential generated by a charge that resides in the 
field of the wave. In the PFR system we have the above cite argument that the 
metric perturbations of the electromagnetic potentials can be ignored. Thus, for 
the description of an electromagnetically bound solid body, we have to take only 
the Coulomb forces into account, and can ignore contributions from magnetic fields 
and electromagnetic radiation. In the PFR system the coordinate distance agrees 
with the metric distance up to order 0(h). Therefore we can state in an invariant 
manner, that the Coulomb potential depends only on the metric distance of the 
particles. We now show how this translates into the TT system. 

We assume that the wave field is slowly varying and the velocity of the charge 
relative to the source of the fields is so small that magnetic fields can be ignored. 
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The corresponding equation to solve for the electromagnetic potential generated by 
a charge q s at rest at x — y = z = is 



d^V~gg^d„A* = ^5\x,y lZ ). (12) 
Our considerations in the PFR system give rise to the following ansatz: 

A° s (r,ct-z) = ^-, r=(x i ) W3 , r 2 = -x^y (ct - z)x* . (13) 
where the Euclidean coordinate distance from the source is replaced by the metric 



distance. We verify that equation (12) is satisfied to order 0(h). For simplicity 



only we consider a wave with the +-mode only where 

r 2 - x 2 (l - /+ (ct - z)) + y 2 (l + f + (ct - z)) + z 2 . (14) 
Using y/—g =1 + (ft 2 ), leaving out terms of order h 2 and higher, we obtain : 



d .P^ha^B - - ( c~ 2 d 2 - d 2 - d 2 - d 2 I - 

x 2 - y 2 \ x y ( z x 2 - y 2 



v r " i 3 ^ - y2 ) 2 ( f > ) 2 1 x2 y 2 f » x2 - y 1 f » , f > d t 

v r 3 + 4r 5 \J+) + 2r 3 + 2r 3 J + + J +° z 



Now the first of these two remaining terms is of order O (ft 2 ) and can thus be 
ignored, the second is proportional to the spatial derivative of the wave field, which 
can be ignored for detectors that are small against the wave length p|. Thus the 
ansatz ( p^ ) solves equation ( |l2] ) except for terms of order O (h 2 ) or proportional to 
d z h. In the same sense the Lorentz gauge holds. The result thus agrees with that 
in the PFR system: The Coulomb potential depends only on the metric distance. 

In the following we assume that this principle may also by extended to other, 
phenomenological potentials, because the time scale set by the gravitational waves is 
by far larger than that of the induced changes of all other fundamental interactions. 
So in the TT system the Coulomb force between two charged particles varies in 
phase with the gravitational wave. It is not hard to see that the local energy 
density of the electromagnetic field, though time-dependent, is only relocated, so 
that the integrated energy density does not change to first order in h, implying that 
radiation effects are at most of second order in h, in agreement with the situation in 
the PFR system. On the other hand accelerated, moving charges produce magnetic 
fields and electromagnetic radiation, but both effects can be ignored for the analysis 
of weak mode excitation in a solid body. 
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5 Many particles 

The main point of this work is to show that a detailed microscopic, many-particle 
model of the detector exhibits two aspects that elude the spring model and the 
normal mode analysis. The first aspect is that the input of energy from the wave to 
the detector occurs only at the surface of the detector. This effect can already be 
modelled with a linear chain , as done in the next subsection. The second aspect 
is the dependence on the nature of the fundamental forces that stabilize a body: 
electromagnetic and gravitational binding forces lead to different mode structure of 
the responde, as shown subsequently. 



5.1 Linear chain model 

In the PFR system TV classical particles with coordinates f^"' = (x, y,z) ,a = 
1, iV obey the equations of motion 

= Y F s , s + ■ r< s > (16) 

dt 2 ^ ss \ ' J 2 dt 2 V ' 

where for simplicity the dot denotes the multiplication of the 3x3 matrix hij with 
the vector f^ s \ and F s > s represents the fundamental forces between two particles. 
Let us consider the special model of a linear chain in i-direction with N = 2n equal 
particles numbered by s' , s — ~n...n, coupled by nearest-neighbor forces 



= -u% (i( s ') -la) for|*'-*| = l (17) 



and a wave with +- mode only, so that the equations of motion are 

Cx^ = -k y ( x (s)- x (s')±i ) + iff± x u (is) 

dt 2 R ^ V V 2 dt 2 v ' 

where k = uj^/m. Without the tidal acceleration, the equilibrium positions are 
= sIq, relative to the center of mass that is identical with the origin of the 
coordinate system. We are now interested in the deviation from the tidal motion 
(??) because the tidal motion itself will be in general too small to be observed itself. 
We set 

e (s) = - ( i + i±] (19) 
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leading to 



f } = -k y ^w_ e (0-(jf«-^))(i+^± fe ) + i^h(fw-iw) 

s'= s±l 

= - k E (& ] -^'^lo^-) (20) 

s'=s±l V 1 ' 

where we have omitted the term \ ^-^r [x^ — X^) because we consider deviations 
from the tidal motion that are proportional to /+ itself, as induced by the wave, 
and thus this term is of order O (h 2 ). Now the following happens: the terms Iq 4^- 
in ( ^o|) cancel for all particles that have two neighbors, only for the particles at 
the ends of the chain they do not. The reason is that under the tidal motion the 
distances between pair of neighboring particles remains constant, thus the induced 
additional forces from the left and right cancel; the pattern of the tidal acceleration 
comes very close to a null mode. So we can rewrite 

t ] =-k y (£ (s) -£ (s) )-W ^(5 s ,„-<Vn). (21) 
s'=s±l 

These equations have a clear interpretation: The deviation from the tidal motion 
is driven by an effective force that applies to the ends of the chain only. Or, we can 
state that it suffices to apply additional forces ^fkl f + /2 to the ends of the chain 
in order to suppress any deviations from the tidal motion. In this case the tension 
of the chain varies uniformly along with the wave strength, but no work is done at 
all against the internal coupling forces. 
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Regarding the energy we have to be careful to use standard expressions, because 
there are metric contributions, as we have seen above, but we can certainly state 
that apart from the tidal motion energy is transferred to the chain only locally at 
the ends. For resonance cross terms in the kinetic energy between the tidal motion 

and the deviations from it play no role, so that ^2 m£ [2 describes the kinetic 
energy of interest. Nevertheless this local picture does not contradict the normal- 
mode analysis in any way. We still have the possibility to decompose the perturbing 
force into normal modes and derive equations for the driving of these modes in the 
standard way. The fundamental mode of the chain is driven strongest, but all other 
symmetric modes are also excited. It is the superposition of all these modes that 
gives rise to the local excitation of the chain. Only refining the spring model to a 
linear chain model could exhibit this property. 

Because we see from (|2l]) that the effective force depends on the microscopic 
property of the chain in form of the equilibrium distance lo, it is necessary to analyze 
a detector on the basis of a more realistic microscopic model that yield information 
on the strength of the driving forces. It turns out that this is also necessary because 
the difference between gravitational and electromagnetic coupling can only then be 
shown. 



5.2 Microscopic forces in the PFR system 

For a more realistic model we can use eq. (|l^) with the fundamental Coulomb, 
gravitational, and repulsive short-range forces inserted. We assume that all these 
forces depend only on the difference vectors — f^°^ and stable equilibrium po- 
sitions -R( Q ) exists. Again we look at the deviation from the tidal motion induced 
by these forces. The transformation 



now leads, assuming small deviations from equilibrium, 

jf«) = £(«) +0(h), 



(22) 



(23) 



to 



.7?(°0 



di 2 ' 
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+ |(h.V)F,,(^ 



(24) 



We will see that this equation for the deviation is identical with that derived for 
the motion in the TT system. The reason is that the transformation (??) basically 
agrees with the coordinate transformation from the PFR to the TT system. For 
the distance we have 



P 



1, 

Oik + 2' life 



-n TT ->j 

-P9ij P 



Oh' 



and thus all potential forces acting on the deviation of the tidal motion in the PFR 
system are identical with that in the TT system. 



5.3 Hamiltonian analysis in the TT system 

We make a potential approximation for the many-particle Hamiltonian, since in a 
fully relativistic approach we had to include necessarily the dynamics of the electro- 
magnetic field in order to preserve energy-momentum conservation. Thus we write 
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the total Hamiltonian for many particles with coordinates 
as 

H = J2 H(S) (25) 

s 

where 

= —p^ s) g ij (ct - zW) p« + V -F ss/ (r< s ) - f< s '\ d - (26) 

2mi V / J z — ' 2 

is the contribution from a single particle. V ss i is the total potential generated by 
the particle s', acting on particle s. As it should be, to each particle only half the 
potential energy is attributed. In the electromagnetic case, the other half, as well as 
the infinite self-energy is subtracted with the contribution from the electromagnetic 
field energy |(| . We assume that the potential is a central potential depending only 
on the metric distance: 

V ss ,(r*^ - r* s '\ct - z< a/ >) = V? sl (yJ- (fW - 9ij (ct - z(*')) (fW - r<«')) J '), 



as derived for the Coulomb potential. 

The evolution parameter is the time t, common to all particles. The Hamiltonian 
(25) conserves the difference between the total energy and the center of mass z- 
momentum: 



= 0. (27) 



But since we have less conservation laws than coordinates, energy and momentum 
transfer from the wave to the particle system has become possible. The change of 
the total energy is calculated using the equations of motion 



2m s 
1 d 



2dz( s 



. s/^s s'^s 



(28) 



that lead us to 



d_ 

dt 



5>w = -Ei^ s)ft 1 d - (s) )^ 

S 8 

\ d M^ s) ~ r* s '\ct-z^) (29) 



2 , 

s,s f s 



where h l i' denotes the derivative, foVsis the partial derivative with respect to the 
second argument only. The derivatives of V s r s with respect to the first argument 
cancel in the sum because of the dependence on f^^ — > only. To first order in 
the perturbation we can approximate 



a ™6 



d 2 V s/s (^ - &'\ct - z^) ~ -V°,' s (r ssf ) ^fK b (ct - z^) (30) 
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where we use i° s , = (r^ — r( s ')) a for short, and r ss , is the unperturbed Euclidean 
distance. Thus we have (with h ab = —h a b in lowest order) 



— — n (s) n (s) - V y 0/ (V ,1 



(31) 



For a pulse of duration T that interacts with the particle system the change of the 
center of mass z-momentum is given by 

AP Z = J2 Kb (d - z (s) ) q$(t)dt (32) 

s to 

(s) 

where q ab is the microscopic contribution from a single particle: 

$(*) = - £ ^ (^) ( 33 ) 

This contribution is quadrupole-like, but it is weighted with derivative of the poten- 
tial. Then the sum or the integral over all particles leads in the case of the Coulomb 
potential to alternating sums where contributions cancel, similar to the summation 
leading to the Madelung constant. On the other hand, if the potential is purely 
attractive, as for gravitation, this effect does not occur, giving rise to a completely 
different response, as well will see. 

For a small detector we can assume that the wave field is constant over the 
particle system represented by the center-of-mass coordinate z(t), so we can change 
the integration variable to r = t — z(t)/c and obtain 

/+oo (a) i.i \\ p + oo 

_ x me t^Mt c - L dThU (CT) Q » i(T) (34) 



with 



Qa b {r)-^ l _ mT))/c . (35) 

For a wave pulse or wave train with duration T, such that h a b (0) = h a b (cT) = 
we use partial integration to write ( |34| ) as 

/+OC 
dTKb (cr) Q' ab (T) (36) 
-oc 

Thus the response of the detector to a gravitational wave pulse travelling in z- 
direction is described by the time-dependence of the microscopic function Q a b- Al- 
ternatively, we may use Fourier transform to arrive at 



AE= / dvh* ab {v)wQ ab {v) (37) 

J — OO 

with the spectral decompositions of the wave: 

h ab (v) = [ h ab (cT)e- 2 ™ T dT, (38) 



and the particle system: 



i(((T))/ C 

Thus Qab represents the effective cross section of the particle system on the micro- 
scopic level. 
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6 Force distribution in a massive body 



We first consider electromagnetically coupled bodies, as ion crystals and metals, 
where the gravitational forces between the constituents play no role. Ion crystals 
have the advantage that all charges can be considered to be point-like and located on 
lattice points. Thus only discrete sums have to be evaluated. As we have seen, the 
forces driving the deviations from the tidal motion in the PFR system arc identical 
to the forces driving the whole motion in the TT system as long as we consider only 
motions of order O (h). 

6.1 Ion crystal 

We first consider a gravitational wave propagating in z-direction incident on a lattice 
of ions. In order for the system to possess a stable ground state, we have to include 
not only the Coulomb potential, but also some (short-ranged) repulsive potential. 
The Coulomb force exerted by particle sf on particle s is given, to first order in h, 

by 

Fc, .* = ^f(f L - ^ + lh J l 'f' X °A , (40) 

to \ ' ssi 1 ssr z ' ssi J 

for the x- and y- directions, the force in z-direction additionally involves the deriva- 
tive of h, which is not of interest here. In the case of the Born-Meyer potential the 
repulsive potential is of exponential form 

V BM (r ssl ) = A aal e-^'" (41) 

with couplings A ss i that depend on the charges. Though it is necessary to include 
this potential, its precise form is not crucial for our further considerations. We 
obtain an additional force from this potential, up to first order in h given by 



( 2 + (3 I hb c x ss/ x ss/ x ss/ 

(42) 

The total force on ion s is given by 

I^s ^ 4TC ° TsS ' TSS ' ' 

~ hab 1^ ~A 3~ + Ass'P X ssl 

+hbc £ ^ + ^ ss,/?£ ^ + ^)) x '" <s ' xls ' 

(43) 

Assuming that the ion chain is in its unperturbed equilibrium state, both the first 
and second terms vanish because the vanishing of the first term defines the equi- 
librium in absence of a gravitational wave, and the second term is just the first 



r BM,ssi 



oVbm, 

Q x a (s) ' 
AssfP 



e -l3r 33 , ! e -/3r 



x<1 S st + A ss iP 
Tss/ ^ 1*ssf 



~A ss /(3 



Tssl 



hab %s 
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multiplied by the matrix h^. Thus the gravitational wave gives rise to a pertur- 
bation of the equilibrium state induce by the force 



AFf = 



^ ( 3q s q s , 1,1. R e-^' ( 1 , p 
^ \ 87T£n r°„ 2 r.., 



sr=£s 



x b 



(44) 



Once the system is driven out of the equilibrium state, we still can ignore the second 
term in (|43|), because for deviations Aa;" s/ ~ h from equilibrium this term is only 
of order h 2 . The first term then describes phonon-graviton interaction. 

The sum over all other ion in ( |44| ) now depends on the dimension of the lattice. 
The sum over the short-ranged part converges even for an infinite lattice, so that 
its contribution is always limited. The sum over the Coulomb part is an alternating 
sum. In one dimension, this sum is always bounded by the first term that is not 
canceled by a contribution from a symmetric neighbor, thus the force has a maxi- 
mum on the endpoints and decreases proportional to i?~ 2 with the distance R from 
the endpoints. In two dimensions we find that the force distribution exhibits the 
correct quadrupole structure |^|]. Therefore acoustic modes are excited. This result 
is obtained only when the short-range potential is included. The Coulomb forces 
alone gives rise to forces that additionally alternate in direction from ion to ion, so 
that we would have arrived at the wrong conclusion that optical modes are excited. 

The forces decay rapidly away from the boundary and do not follow the linear 
law of the tidal accelerations. Thus, the deviations from the tidal motion excite a 
broad spectrum of modes, of which few are resonantly driven. We now present a 
general argument that this pertains to the relevant case of three dimensions. 



6.2 General case 

Since the temperature of the body must be low in order to achieve the desired 
sensitivity of a detector, we assume that the mass is a perfect crystal, either an 
ion crystal with discrete charges locate on some lattice, or a metal with ions on 
some lattice and the electron gas in between. The body is decomposed into a finite 
number of elementary cells Vi with their charge centers at fi , that are (i) electrically 
neutral, 

/ d D rp(f) = (45) 

JVi 

and (ii) do not possess an electric dipole moment, 

/ d D r (f- n)p{f) = 0. (46) 

JVi 

D is the spatial dimension of the lattice. We now consider a charge element q' = 
p{f)dV located at f and the perturbational force exerted on it by the elementary 
cell V located at fi = 0, according to (ff"I|). For If] ^S> 1//3 we can ignore the 
short-range potential and have 

ma if jD 3/j(r) (r — f) . ^ ,b _»c f .-s 

F s (r)=q / d r — — — -g-^ r - r) (r - r) . 47 

Jv 87r£ o IF - rf 

We expand the integrand into powers of f = ( x, y, z) up to second order. Then 
the integrals of the first and second order vanish due to conditions ([l5]) and fl46|), 
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respectively. As an example, if h is of +-polarization, the forces in two dimensions 
are given by 



where the integrals 



x' (6x' 2 - V 2 ) xl (-3x' 2 + 12y' 2 ) y' (l2x r2 - 3y' 2 ) 
2^7 h + 2^7 h + ~i /3 



h 



,d 2 3p(r) 
a r x 



87T£o 



d D ry 2 



Mr) 
8irso ' 



r xy 



87T£ 



F, 



< 



E 



Q>,Vi 



< const 



E 



i 



(48) 



y' (I2x' 2 - 3y' 2 ) y' (-9a;' 2 + 6y' 2 ) x' (-3x' 2 + 12y' 2 ) 
2V7 h + 27* h + ?f h 



(49) 



(50) 



describe the quadrupole moments of the electric charge distribution in the elemen- 
tary cell. The generalization to three dimensions is straightforward, with the same 
qualitative properties: The force exerted by some elementary cell on a charge el- 
ement q' always decreases at least proportional to r' -4 with the distance between 
the charge and the cell. If the quadrupole moments vanish, as for cubic lattices, 
the decrease is even faster by two powers of r'.This implies that for a given charge 
element q' the sum over all elementary cells always converges in D — 1,2, or 3 
dimensions. Therefore the force on any ion in the body is bounded independently 
from the size of the body, 



(51) 



Further, the reflection symmetry of the lattice implies that all forces from cells in 
a volume that possesses reflection symmetry around the charge element add up to 
zero. Hence also in the general case the charge element feels a force that depends 
on its distance R to the surface of the body. A crude estimate gives 



< 



E 

,\rl— Pi \>R 



F y 

q',Vi 



(r - n) 



< const 



1 



R 



D-4 



(52) 



\Pi-Pi\>R 



Depending on the dimension, we observe that the force decreases with the distance 
from the surface, at least with l/R in three dimensions, if the quadrupole moments 
Ii,h,h all vanish the decrease is even stronger. Scaling with the lattice constant 
a leads us to 



< 



F 



q> 



(I) 



i-D 



(53) 



Because the sum over the forces from the short range potential decreases even faster, 
the total perturbational forces, which are maximal on the surface, decrease to about 
10~ 3 of the surface value within 1000 atomic layer, which is about lmicrometer. 
When we integrate the forces (|^,|^|) over an elementary cell in order to obtain the 
mean force on the cell, we again loose two powers due to ( fl5|) and (|4^), resulting 
in mean forces between two cells Vi and Vj that decrease at best proportional to 
\rj — rip 6 . Thus we conclude that the bulk of the material remains, apart from 
tidal motion, unaffected by the forces induced by the gravitational wave. This 
result has its origin in the nature of the electromagnetic coupling with its charges 
of different signs, and the structure of solid bodies with reflection symmetry of the 
elementary cells. 
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6.3 Gravitationally coupled matter 

Clearly, the result of the preceding section was due to the conditions (45 4^), but 



a G (?) = I d D r "gl ±—^h bc (f - ff {? - f) c , (54) 



(45) does not hold for the attractive gravitational forces. The acceleration of a test 
mass is given by 

D„ P, n(r) (g - r) 
G |/* — f\" 

where p m is the mass density of the body. Assuming that the wavelength of the 
gravitational wave is large compared to the dimensions of a homogeneous massive 
body, we can take p m and h to be constant and are able to evaluate the integral for 
simple geometries. For a +- polarized wave propagating in z-direction the maximal 
acceleration on the surface is calculated to 

p m 7rr 1 . 1 | for a cylinder of radius r and lenght I 

^"^ (l+tan 2 1) 3/2 ^ ^ G ^ P °^ a ° 0ne °^ °P enm S an gl e 4> an d height I 
Pm^if at the surface a sphere of radius R. 

(55) 

Naturally, there exists an angular dependence of the forces of quadrupole charac- 
teristic. Thus the forces grow linearly with the linear dimensions of the massive 
body. The maximal acceleration on the equator of a sphere is given by 

\MlL = \\h\9 (56) 

where g is the surface acceleration of the mass. This force will truly excite the 
quadrupole modes in the bulk of a body and is able to do work against the grav- 
itational and electromagnetic forces that keep the body together. Comparing this 
with the tidal acceleration in the PFR system, 

\^C F J=^ 2 R\h\, (57) 

we see that ([36|) is several orders of magnitude smaller than (|57f). In the PFR 
system the forces drive primarily the tidal motion, only accelerations of order ( [56l) 
drive the deviations. 



7 Detector Types 



7.1 Rotational Detectors 

From the structure of the microscopic quantity Q Q b(r) we can immediately identify 
the different types of detectors. If the body is rotating relative to the TT frame 
with some frequency u, then the momentum contribution 

T,^P« ] Pb S) ~" 2 co S 2a,i (58) 

s s 

dominates and gives rise to a response of order h. This type of detector was first 
suggested by Braginsky ^ ||. We estimate the response of a rotating mass to the 
gravitational wave using our model. From (]34|) we obtain a momentum input to the 
center of mass of 

dP z /dt — h cos (coot + ip) cos2o;t tJ — (59) 

24c 
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where I is the length of the bar, M its mass, and ho, luo, and tp are the amplitude, 
frequency and phase difference of the gravitational wave, respectively. This is for 
the case where the wave vector is perpendicular to the plane of rotation. In the 
ideal case lo ~ ujq/2, tp = or n, the mean energy input is given by 

E = ±h ^ r . (60) 

Note that the change can be of either sign, thus the gravitational wave cannot only 
be absorbed, but can also stimulate the emission of gravitational waves from the 
system. For reasonable values (bar of lm, 1 — lOOfcg mass, oj/2it ^ 10— 1000 Hertz, 
ho ~ 10~ 20 ) the attainable energy input ranges in about 

E ~ Kn 13 ...lCr 9 W. (61) 

Though this is quite large compared to the resonant detector, it seems questionable 
whether this change in the rotational energy of the bar can be measured. Certainly 
the acceleration of the rotating bar in the direction of the wave, given by (|59|), is 
too small to be detectable. This gives additional justification to the negligence of 
the acceleration of the center of mass in the transformation to the PFR system Q . 

7.2 Resonant Detectors 

For a non-rotating mass, the lowest order contribution stems from the time deriva- 
tives of the positions and momenta induced by the wave and thus leads to a response 
of order h 2 . This is the Weber detector. Our considerations have shown that the 
energy input into such a detector occurs primarily at its surface. This implies that 
when resonance is discussed, we have to take the time into account that is needed 
to transport form the surface to the interior of the body. For a material with high 
Q-factor the velocity of transport is given by the velocity of sound. Thus it takes a 
time 

T v = A (62) 

where L is the diameter of the body and v s the velocity of sound before energy 
reaches the center of the body. This time limits the onset of resonance. For example, 
for GRAIL [?] with diameter L = 3m, v s ~ 4000to/s we have T v « 3/8 milliseconds. 
Thus this time scale will play a role for the detection of millisecond pulses. In order 
to arrive at the resonant amplitude, pulses must be considerably longer than T v . 
How many oscillations are needed to arrive at the maximal resonant amplitude, can 
also be estimated in the normal-mode model. 

In general, if a weakly damped oscillator characterized by a frequency u>o is 
excited by a force F — euj 2 sin (urb) near resonance, uj rs u>o, then if the oscillator 
is initially at rest, the subsequent maxima of the amplitude of the oscillator will at 
the beginning follow a linear law 

2 

\A max \ ~ ^-t (63) 
2w 

until the resonant amplitude of order C — eCo is reached after a transient time 

2Cq 2Cq 

Ttrans ~ ~ ■ (64) 

LO UJo 

Thus e.g. it takes for a signal of frequency u> = 1000Hz about 2s until it is enhanced 
by a factor Cq — 1000. Thus 2000 oscillations are needed before the maximal 
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resonant level is reached. In general the effective cross section of a resonant mass 
detector is proportional to (Co) 2 , but for pulses shorter than T tra ns an additional 
factor (T /Ttrans) must be taken into account. This is relevant to the detection 
of millisecond pulses with a detector like GRAIL that operates at a fundamental 
frequency of 650Hz . 

A further point to consider certainly is the impurity of the material. When 
grains of material stick together, we also have to consider internal boundary that 
behave like surfaces. But because this results in random effects, the impurity will 
hardly improve the behavior of the detector. 



7.3 Microscopic Detector 

We have shown that in the TT system the metric effects on the electromagnetic 
coupling cancel in the bulk of a massive body, in the PFR system the deviations 
from the tidal motion are driven only at the surface. This raises the question 
whether this type of surface effect could be used in some other way for detection 
of gravitational waves. Similar to expression ([56]) for the gravitationally coupled 
matter, we can estimate the acceleration induce by the gravitational wave on the 
surface of an ion lattice by 



\ a EM\ mav 



Z + Z- 



m±d 2 Trso 



(65) 



where Z +1 Z- and m + , m_ are the charges and the masses of the ions, respectively, 
and d the lattice constant. A geometric factor of order 1 has to be included in 
addition. This factor will depend on the structure of the lattice and the precise 
behavior of the repulsive potential. The numerical value of ( |65| ) is of the order of 
h ■ 10 16 ms~ 2 (for KC1) which is several orders of magnitude larger than (|57j). This 
suggests that the piezoelectric effect might be used for experimental verification, in 
a similar way like a pressure sensor works. For h in the order of 10 -22 accelerations 
are of order 10 -6 ms -2 and possibly are within the reach of sensitive detectors. If we 
consider a piezo crystal under strain our analysis has to be revised. Because we then 
have a dipole moment in each elementary cell, the leading order of the forces will 
be proportional to r'~ 3 as compared to r' -4 in ( fl8| ) and (^). Due to the uniform 
direction of the dipoles and the broken refection symmetry the contributions in the 
integration will no longer cancel, so we expect that the forces grow logarithmically 
with the dimension of the crystal and the excitation occurs truly throughout the 
bulk. We hope to be able present a detailed analysis of this case, together with 
possible detector design, soon elsewhere. 



8 Conclusions 

We analyzed the resonant mass gravitational wave detector form a microscopic point 
of view, using the wave guide (TT) frame of reference along with the conventional 
PFR system. In the TT system the variation of the Coulomb field of the constituent 
charges of the body gives the dominant contribution, the resulting forces agrees 
with those driving the deviation from the tidal motion in the PFR system. For an 
electromagnetically coupled body with reflection symmetry, this force distribution is 
such that the forces are restricted to a small surface layer, with no bulk force. Hence 
the relevant energy input occurs at the surface only. This effect has its origin in the 
fact that the pattern of the tidal forces comes close to a null mode of the system, 
not doing work against the coupling forces, as we saw from a linear chain model. 
This result is not in contradiction to the standard normal mode analysis, rather 
it reflects a local property of the response compared to global nature of normal 
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modes. For the gravitationally coupled body, we observe a linear force law and 
bulk excitation, but the part that causes resonant response is orders of magnitude 
smaller than expected from the tidal forces in the PFR system. 

Regarding resonant mass detectors, the local nature of the driving force makes 
it necessary to consider the time scale on which the energy is transported into the 
interior. This time scale is set by the velocity of sound. As a result, it might take too 
long a time for a resonant detector to be excited measurably by short millisecond 
pulses as are emitted by collapsing stellar objects. 

Finally, we presented ideas for a new type of microscopic detector that employs 
the piczo effect. The force pattern on a crystal is basically that of a quadrupolc- 
like distributed pressure change, the surface force is independent of the size of the 
crystal. For a crystal under strain the force even acts throughout the bulk. An 
analysis how the desired sensitivity can be achieved was outside the scope of the 
work presented here. 
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